setwd("~/Downloads/dataverse_files-1")
library(ergm)
library(stargazer)
library(coda)
library(latticeExtra)

sink("PSRM_selection_model_log.txt")

# Selection model: table 7 
# Making the selection bias table
selection = read.csv("PSRM_selection_bias.csv", sep = "\t")
mod <- glm(responded~., family = binomial(link="logit"), data=selection)
stargazer(mod, title = "Selection Bias Model", out="results_from_table_7.txt")


sink()
